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Abstract 

We continue our previous work on the flavour-conserving leptonic decays of 
the Z boson with neutral heavy leptons (NHL's) in the loops by consider- 
ing box, vertex, and self-energy diagrams for the muon decay. By inclusion 
of these loops (they contribute to the input parameter Myy) we can probe 
the full parameter space spanned by the so-called flavour conserving mix- 
ing parameters ee m i x , fJ-fJ-mix,TT m i x . We show that only two diagrams from 
each class (box, vertex and self-energy) are important; further, after renor- 
malization only two box diagrams 'survive' as dominant. We compare the 
results of our analysis with the existing work in this field and conclude that 
flavour-conserving decays have certain advantages over traditionally consid- 
ered flavour- violating ones. 
PACS number(s): 14.60. St, 12.15.Ff, 13.35.Bv 
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I. INTRODUCTION 



We have previously considered a simple extension of the standard model (SM) with 
an enriched neutral fermion spectrum consisting of a massless neutrino and a Dirac neutral 
heavy lepton (NHL) associated with each generation [@-^]. Several parameters can be used to 
characterize the model: 'flavour-conserving' mixing parameters ee mix , fifi mix ,TT mix ; 'flavour- 
violating' mixing parameters e/i m j X , er mix , \iT m i X and the mass scale of NHL's (assuming 
three degenerate NHL's). We considered the effect, via these parameters, of NHL's on 
flavour- conserving Z boson decays to charged leptons and on the W boson mass, My/. 
However, in our earlier work, we neglected all mixing parameters except rr m j X , which is the 
least well constrained. Here we generalize our analysis by considering the case of arbitrary 
mixings ee mix , ii\i m ix and TT mix . Our previous neglect of ee mix and \x\x m ix allowed us to also 
neglect a number of contributions to the muon decay corrections which feed into Mw as an 
input parameter. Including these couplings, non-SM box, vertex and self-energy diagrams 
contributing to the muon decay (see Figs. [I], ^, [5D may become important for the calculation 
of Mw- In our previous paper as a result of the assumption ee m i X = nn m ix = 0, 
only oblique corrections (corrections to the W propagator) had to be considered. Here 
we consider the full set of corrections. Still, we assume here vanishing flavour-violating 
mixing parameters: efi mix , er mix , fiT mix = 0. These parameters, if nonzero, lead to further 
complications, which in general require, as argued in a recent work ||, the renormalization 
of the mixing matrix. This is an interesting topic by itself; nevertheless, it is not crucial for 
our considerations. We note that the smallness of e[i m ix is confirmed by experiment 

The inclusion of the arbitrary flavour-conserving mixing parameters completes our stud- 
ies of the NHL's impact on the processes considered here. We compare our constraints on the 
parameters of the model with those coming from the traditionally favoured flavour- violating 
processes, such as /i — > ej, r — > eee, Z — > e/x, etc P,|7Hl4f. We find that the processes we 
consider have certain advantages over the latter ones. 

This paper is organized as follows. In Sec. II below, we briefly review a superstring- 



inspired SU(2)l x U(1)y model of neutrino mass and the constraints on the mixings and 
masses of the model. In Sec. Ill, we present the additional muon decay corrections, iden- 
tifying which contributions are important. Ultimately, our earlier results can primarily be 
improved by the tree level modification of the vertex by mixing factors. In the limit of large 
M^, only two box diagrams finally contribute but these are numerically only marginally 
important. Given the muon decay corrections, we also present the one- loop modification of 
the constraint on TT mix . In Sec. IV, we consider more generally the work done in this field. 
We contrast the sensitivity to the presence of NHL's in flavour- violating processes with the 
results for flavour- conserving processes. We include a calculation of the flavour-violating 
leptonic decays of the Z boson in our model. We summarize in Sec. V. 

II. A SUPERSTRING-INSPIRED SU(2) L x U(l) Y MODEL OF NEUTRINO MASS 



Here we briefly describe the model of neutrino mass which we consider. For more details, 
we refer the reader to the original papers or our previous work . The model extends 
the neutral fermion sector of the SM by two new weak isosinglet neutrino fields Sl) per 
generation. With total lepton number conservation imposed, the mass matrix is given by 
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Each ul, tir, Sl represents a collection of three fields, one for each family. D and M are 3x3 
matrices. The diagonalization of the mass matrix yields three massless neutrinos (vi) along 
with three Dirac NHL's {N a ) of mass Mn ~ M. The weak interaction eigenstates (z//,Z = 
e, n, r) are related to the six mass eigenstates via a 3 x 6 mixing matrix K = (K^, Ku): 



">= £ (K L ) H u tL + £ (K H ) la N aL . 

i=l,2,3 



(2) 



a=4,5,6 

The mixing factor which typically governs flavour- conserving processes, ll m i x , is given by 
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nnx = £ ( K H)l a {K ] H ) al ; 
a=4,5,6 



I = e, fI,T 



(3) 



and the flavour- violating mixing factor 11 mix is defined as 

it mix = E ( K H)ia(K ] H ) al > ; l,t = e,(l,T, (4) 

a=4,5,6 

Further, the following important inequality holds 

\H mix\ ^ Umix I I mixi I ^ I . (5) 



This implies that one might observe nonstandard effects in flavour-conserving processes even 
if they are absent in flavour-violating processes. 

We note here existing constraints on the parameter space of the model. Indirect con- 
straints on the flavour- conserving mixing parameters ee m i X , fifimix, TTmix have been obtained 
from a global analysis of results including lepton universality measurements, CKM matrix 
unitarity tests, W mass measurement, and LEP I measurements. These constraints arise 
primarily at tree level due to the modification of couplings from those of the SM. Nardi et 
al |15| have found the following upper limits: 

ee mix < 0.0071 

W m ix < 0.0014 (6) 

TTmix < 0.033 

Since the limit on the parameter rr mix plays (as the least stringent one so far) the most 
important role in our analysis, we will pay further attention to its source. The /i — r 
universality test is based on the r leptonic decays compared to the \i leptonic decays, with 
the result given as the ratio of the couplings of r and /i to the W boson, g T /g^. The tree 
level ratio is found from 

r(r -> evv)/T SM {T -> evv) _ _ 1 - TT mix 

This measurement has undergone substantial improvement over the last while. With the 
most recent result g T /g^ = 0.9994 ±0.0028 [I6|, the constraint on TT m i X is improved from its 



previous value of 0.033 by a factor of about three. To reflect this improvement we present 



most of our results either for the values of TT mix ranging from 0.033 to 0.01, or in a general 
form with TT mix as a variable. In a few cases (e.g. when quoting results of others on flavour- 
violating processes), however, we only use TT mix = 0.033. Finally, we note that these indirect 
limits depend very weakly on M^', this point will be illustrated at the end of Sec. III. 

Since NHL's have not been directly observed in the Z decay Z — > Nu, we focus on 
NHL masses 3> Mz, Myy, Mfj- These can be probed indirectly via loop effects in either 
flavour- violating or flavour-conserving processes. As argued in our previous work |]J , only in 
this case are the contributions of NHL's via loops possibly significant, due to the violation 
of the Appelquist-Carrazzone decoupling theorem |17]]- Analogous to the behaviour of the 
top quark loop contributions in the SM, quadratic nondecoupling (amplitudes ~ M^) often 
results here. 



III. NHL CONTRIBUTIONS TO MUON DECAY 
A. Box diagrams 

We first consider the box diagrams contributing to the muon decay, as depicted in Fig. |I[ 
Diagrams of Fig. |I|b,f,g each come in two varieties, with either massless neutrinos or NHL's 
in the loop. All diagrams without NHL's are similar to their SM counterparts; the only 
slight difference comes from the mixing factors in vertices (such as 1 — ll m i X , see Appendix 
A). All the box graphs are finite. 

The results of our computation of the diagrams of Fig. |l]a-g are given in Appendix A. (The 
QED box amplitude of Fig. Pi, M. ie w^i is given in Ref. ||18|| .) The dominant nonstandard 
contribution in the limit of ^> Mz, Mw, Mh comes from just two graphs depicted in 
Fig. |l|e, one with Higgs H and one with neutral unphysical Higgs x- To show that these two 
graphs are dominant we could take the large limit of the exact results given in Appendix 
A. We would find that only these two graphs exhibit quadratic nondecoupling, i.e., quadratic 
overall dependence on Mjy. The remaining graphs with NHL's are either constant in the 



large limit, or decouple as 1/Mjy-. However, here we prefer a more intuitive approach 
based on dimensional analysis considerations and power counting. 

The amplitude for the diagram with the Higgs boson H (Fig. |l|e) is given by (we sum 
over NHL's N a , Nb with M^ a = M^ b = and neglect external momenta in the internal 
propagators) 

ka r dk 4 _ -ig 2 M N , 1 + 75 i ig 2 M N * 
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Various mixing factors can be collected as 

k mix = {K{K H ) ia {K ] H ) aii {K H ) eb {K ] H K L ) b0 = (Kl) ilx (K L ) ej ee mix ^ mix , (9) 

where we used efi mix = fir mix = 0. Neglecting some constant factors which we will restore 
later, we get 

M^hn ~ < / {k2 _ M%)2{k2 k _ M ^ )(k2 _ M%) . (10) 

Note that the Lorentz structure of the amplitude is such that NHL propagators y_ l MN 
contribute as k 2^ M -i rather than pz^z ■ In the limit of large we can neglect all masses 
and momenta except Mjy, obtaining 

d 4 k 



M, NHN ~M N j w4 {k2 _ M%)%k2 . (11) 

The integral is expected to be of the form (Mn) p ; power counting yields p = —2, so indeed 
the amplitude depends quadratically on M N : 

M,, x „x - .Uv.\/ V 2 = M 2 N . (12) 

We can further improve our estimate by restoring the constants collected from Eq. [8], 

■M-4>NHN = cMtTee—2~lT7^ ee rriix^^mix-, (13) 
S W M W 



where 

■ 2 

Mtree = "g^EM 1 + ^^[v^ + Js)^) (Ki) ^(K L ) ej , (14) 

is the tree-level amplitude. The remaining numerical factor c can be found from the exact 
result given in Appendix A. It is equal to c = 1/647T. 

Similarly, the amplitude M^NxN is, in the large M N limit, equal to M^nhn- Dimensional 
analysis can also be applied to the remaining boxes, confirming that J^A^nhn and M.^n x n 
are the only box diagrams with quadratic nondecoupling. 



B. Vertex diagrams 

We next consider together vertex corrections and corrections to the external charged 
leptons. Diagrams modifying the W ' [iVi vertex are depicted in Fig. |2|. Another set, one that 
modifies the Weuj vertex, is not shown. 

The sum over the depicted set of diagrams gives the muon vertex amplitude M. vertex'- 

■^vertex = -M-fiuZ + -MfiNZ + -M-ZWfi + ■MyW^ + -M-WZu 

+ M W ZN + MfiZN + M W HN + MjHN + MfaN (15) 

Explicit expressions for each of these amplitudes are given in Appendix A. They are divergent 
and we renormalize them with the SM form counterterms Jl8| (renormalized quantities are 
distinguished by the hat): 

A" = W + SZY-SZ™ + 6Z£, (16) 



where 



«F-«? = -^(f-7 + 1n*-ln^) = I") 



^=-EE(mJ) + ^(21„^-l), (If 



2/e with e — ► is the pole of the dimensionally regularized amplitudes and A is the regularized 
photon mass. = T)^ N + + £^ + M + is the left-handed part of the muon 
self-energy, with the individual terms corresponding to the loops shown in Fig. |^. All these 
contributions are given in |1J and The term which we use specifically below, E^, 

is given in Appendix A. In our scheme, the renormalized charged lepton self energies do 
not contribute directly, but rather through the renormalization constant SZ^. Cancellation 
of divergences occurs as usual between the vertex loops and the counterterm contributions 
but we focus here on the dependent terms only. Looking for the dominant graphs 
in the limit 3> Mw, Mz, Mh, we find (either by taking the limit of exact results or 
using dimensional analysis and power counting) that the graphs of Fig. 0: have quadratic 
nondecoupling. However, both infinite and finite parts of these two graphs are cancelled in 
the large Mjy limit by the term (see Fig. |3|b) in the counterterm 5Z£- Therefore there 
remain no Mjj dependent terms in the renormalized vertex diagrams. 

This curious cancellation can be seen either explicitly (both and M^hn, M-4>xN are 
given in Appendix A) or, better yet, after applying the symmetries of the theory. The way 
to go is to study the more familiar case of a 7// vertex. This vertex is modified from its 
tree-level value iej/j, by the one-loop diagrams shown in Fig. |] as (we show only vector and 
axial vector corrections) 

ie-ff, -> iej^l + F v ) - ie^^ 5 F A . (19) 



We now use a Ward-identity ||18|J19| , which relates the vertex formfactors F Vi a evaluated at 
(Pi +P2) 2 = (pi,P2 are lepton momenta) to charged lepton self-energies represented by the 



counterterms 5Z\ 



V,A- 

where 6Z l v = \{5Z l L + SZ l R ), 5Z l A = \{8Z l L - 5Z l R ), and £ 7 z(0) = § ^A Mff is the term 
originating in the bosonic loops of the 7-Z mixing. At small Mn the graphs with unphysical 
Higgs <f> are negligible, however, with M N ^> M w , M z , M H two types of graphs dominate 



the left-hand side of Eq. |20|: the irreducible vertex (formfactor) Fy^ (see Fig. |]c) and the 
self-energy (its vector or axial- vector part) 5Z V N A (Fig. |3|b). Since the right-hand side of 
Eq. is not affected by the NHL's, it remains constant and (very) small with respect to 
5Z V N A at M N = O(TeV). Hence the only way to meet the above formula is to have 

F^ A N + SZt N A = (21) 

in the limit of large Mat. If we now return from the jll vertex to the W \iv vertex, we find 
a similar result (for proof see Appendix B): 

A^hn + A^n + 6Zf = 0, (22) 



that is, the two dominant nonstandard contributions from Eq. [R] cancel exactly, includ- 
ing the finite parts. Since the remaining nonstandard contributions are not enhanced by 
the quadratic nondecoupling and are suppressed by the mixings, vertices can be reliably 
represented by the SM terms. 



C. Neutrino self-energy and its renormalization 

Half the neutrino self-energy diagrams contributing to muon decay are shown in Fig. [5]. 
The corresponding self-energy is denoted as S"* 4 . The other half consists of the same loops 
sitting on the bottom neutrino leg with the corresponding self-energy E" 8 . In all these 
diagrams, we sum over the internal massless neutrinos Uk,k = 1,2,3. In principle, the 
graphs with v k replaced by N a are also present, however, they are suppressed by the large 
mass Mm- 

The unrenormalized neutrino self-energy YT l (I = e, /x) has the form 

TP = \t% ^(1 - 7b), (23) 

where receives contributions (given in Appendix A) from the diagrams of Fig. [5[ The 
amplitude for those diagrams, in terms of S^' , can be shown to be equal to 



M sd f = -Mtree^i (24) 

where the factor ~ comes from our dealing with the external wave function rather than the 
neutrino propagator. 

Let us now investigate the question of the renormalization of In this case the 

counterterms are modified from their SM form. The problem is how to renormalize a part 
of a theory where interaction eigenstates are different from mass eigenstates. Curiously, this 
also happens in the SM quark sector [p0|| . The difference is that in the SM the problem is 
circumvented by arguing the off-diagonal quark mixings are too small to have any effect in the 
loops and the renormalization procedure is effectively simplified to that of mass eigenstates 
being also flavour eigenstates. In our model, we cannot neglect the 'off-diagonal' mixings 
(their role is assumed by ll m i X ), since they (in combination with TeV NHL masses) lead to 
the dominant terms in the predicted deviation from SM results. This problem was studied 
in Refs. In Ref. || it was shown that in general the renormalization of the divergent 

amplitudes requires the renormalization of the mixing matrix. In our model, the amplitudes 
can be renormalized without the renormalization of the mixing matrix, if the assumption of 
zero flavour- violating mixing parameters is made. Our scheme is a straightforward extension 
of the SM counterterm. 

We start with the counterterm Lagrangian, which has the same form as that of the SM. 

i SZ e L V e jdu e + i SZ^Vf, jdv^ + i 5Z T L V T fiv T . (25) 

Weak eigenstates v\ are given in terms of mass eigenstates Ui, N a in Eq. |[ This gives us, for 
the product 7^z/ z , 

Vm= £ V i {Kl) il {K L ) lk v h + ...{v i N,Nv k ,NN), (26) 

fc,i=l,2,3 

and Eq. ^ thus contributes the following massless neutrino counterterm: 

E {SZl(K[) ie (K L ) ek + SZ^Kl)^)^ + 6Zl(K[) lT (K L ) Tk }v z M- (27) 

fe,i=l,2,3 



In our case we sum over internal v\. but not over external v^. The graphic representation of 
the relevant counterterm (embedded in muon decay) is in Fig. |6]. 
The amplitude for this diagram is 

Mo = -\m™ E E ^i{K\UK L UK{) k ^K L ) ey (28) 

l=e,fi,r fc=l,2,3 

Again, the factor | comes from our dealing with the external wave function rather than 
the internal propagator. Mf^ee ls ^ ne ^ ree l eve l amplitude for muon decay in the SM. The 
mixing factors (K* L ) k and (Kl)^ originate at the fiWuk and eWuj vertices, respectively. 
The amplitude M. c can be further simplified, 

Mc = -\mH{ E SZ l L (Kl) a E (K L ) lk {K{) k ^K L ) ej 

~~ l=e,fi,r fc=l,2,3 

= -\m%Z E SZ l L (Kl) u (5 ¥ -li, mix )(K L ) ej 

l=e,fi,r 

= -\M™6Zt{l - w mix ){K{)^K L ) e3 

= -Uzt{\ - w mix )M tree . (29) 

The factor {K ] L ) v {K L ) ej was absorbed by M tree = M s t ^ e {Kl)^{K L ) er 

Now we can write down the final expressions for the renormalized amplitude Xi se if and 
the renormalized neutrino self-energy : 

Yf l 5Z l 

M S elf = Mself + M C = Y Mtree - U mix)M t ree, (30) 

±2 = ^+SZ l L (l-ll mix ). (31) 

The constant SZ l L was given in Eq. [H| 

To prove the cancellation of the infinities, we note that the infinite part of 5Z l L is given 

by m 



8Z %° = "Si^ + i + z^^K (32) 



where A M = - — 7 + ln47r + ln/z 2 . The infinite part of the neutrino self-energy is (see 
Appendix A) 



E^ a 2 { A /l/f2 ^wia;(l ^mix) "I - r>(l ^mix) a 2 ^w*iit(l H"mix) 

+ -^(l-// m „) 2 }A M . (33) 



From the formulae above it can be easily seen that infinities cancel out in Eq. [3l. 

We now investigate the large behaviour of the renormalized neutrino self-energy E^ , 
this time using exact results. The two diagrams of Fig. |5]c contribute to the self energy 
with an overall factor of M^. For large M^, the coefficients of these diagrams contain the 
functions 

B (p; M H , Z>W , M N ) ~ 1 - 2 lnMjv, 

B^Mh^Mn) ~ -0.25 + In M*. (34) 



This implies quadratic nondecoupling for E|f (p) and E£(p) such that 



2 r 



(p) + ^l(p) — 2n 4s 2 ^ mix ^ Mmix) j^2 



1 . 3 , ' 
-A„ + --lnAr* 



(35) 



E^ (see Fig. ^b), which contributes to E^ via the counterterm (see Eqs. [18], [!l]), is 
given in Appendix A (Eq. |A8| ). From here we can see that (once again) E^ not only cancels 
infinities in E|;(p) and E*(p), but, in the large Mat limit investigated, it also cancels the 
finite parts. As a result, there is no quadratic nondecoupling in the renormalized neutrino 
self-energy and, as in the case of irreducible vertex corrections, it suffices to consider just 
the SM loops. 



D. Results 

The loop corrections to muon decay modify the quantity Ar in the implicit relation 
between My/ and G M as follows |J 

2 2 7ra ( 1 1 \ 

M w s w = ^ _ (1 ~~ 2 e<2mi:c — ^HH-mix), (36) 

where 1 — \ee m i X — \\i[i m ix is the tree-level correction in our model and Ar can be written 

as 



Het w (0) 

Ar = 172 + S v- (37) 

m w 

Svy(O) is the renormalized self-energy of the W boson which we previously calculated 
The parameter 5y is the sum of the boxes, irreducible vertices and self-energies calculated 
in the previous sections, along with the equivalent contributions to the Wev vertex, 

Sv = MjeV^ + Mbo, + fa + fa _ 1^ _ -U^ (3g) 
Adtree 2 2 

Based on the previous sections, we expect that 8y can be reliably represented as 
o v — o v + o b — o v +— — —-rj^ee m i X ^L^L miX) 



where 5y M is the SM value [O] and the rest comes from just two box diagrams (Fig. 

Numerical results for the corrections to muon decay are shown in Table |. As input data 
we used the following set (henceforth the standard set): Mz = 91.1884 GeV, a" 1 = 137.036, 
A = = 37.281 GeV, M H = 200 GeV, m t = 176 GeV. The mixing parameters used are 
ee m j X = 0.0071 and /i/i m j X = 0.0014 while for TT mix we show results for both the minimal 
and the maximal value allowed, and 0.033, respectively. rr m j X is the least well constrained 
mixing however there is no particular theoretical motivation to assume that it is actually 
larger than the other mixings. Hence we give results with rr m j X suppressed in order to 
illustrate the dependence on ee m ix and mi m ix- For rr mix = in the first three lines of the 
table we show the contributions of the self-energy, vertex and box diagrams to 5y (line 4) 
for NHL masses Mn of up to 30 TeV. Also shown (lines 5,6) are Svf(0)/M^/ and Ar since, 
ultimately, we are interested in NHL effects in the observable Mw (line 7). The SM values 
are given in the first column. The results confirm expectations from the previous sections. 
There is no nondecoupling for self-energies and vertices and there is a quadratic dependence 
on Mtv, in the large limit, for the boxes. The boxes are becoming important at very 
high masses. Still, they are small compared to the change in T,w(^)/M w . This is due to 
the fact that the dominant boxes enter with the coefficient ee m i X fifi m i X (see Eq. [13]), while 
the correction to the W propagator is proportional to k HH = e^ mix + ^fi^ix + TT m,ix H]> 



which is allowed to be larger given the current bounds on the mixings. The W mass jumps 
from M^ M = 80.459 GeV to M w = 80.537 GeV at M N = 0.5 TeV mainly as a result of 
the tree-level correction factor (l — \ee m i X — \^^ m ix) (see Eq. [36|). After that it rises very 
slowly until the M/v dependent amplitudes become dominant above 5 TeV. 

The results for rr m j X = 0.033 case (Table |, lines 8-11) are similar. We only show 
5y, £w(0)/M^, Ar and M w since the boxes, self-energies and vertices change slowly with 
TT mix (they only depend implicitly on TT mix , via Sjy), as illustrated by 5y in line 8. The 
relative impact of nonSM boxes (Sy) on Ar compared to that of Svy(0)/M^ decreases with 
increasing TT mix . Sjy(0)/M^ ~ rr^Mff corrections actually violate (for > 5 TeV) the 
perturbative unitarity bound discussed in Ref. [1] . 

To sum it up, the analysis of Ref. [|IJ turns out to be basically valid even after the 
restriction ee m i X = fifi m i X = is relaxed. The numerical predictions can be improved 
by the inclusion of the tree-level correction (l — \ee m i X — ^^mi^j, while the largest loop 
corrections, the box diagrams of Fig. |l|e, are only marginally important. Only in the case of 
the Z decay into e + e~, with ee m i X and TT mix now made comparable, does the character of 
the M N dependence change (see below). 

The impact of these results is shown in Figs. and We give the Z leptonic widths 
as a function of NHL mass for the two cases, ee m ix = and ee m i X = 0.0071, (fifi m i X is 
negligible) in Figs. and [|, respectively. These figures show the widths for three values 
of TT m i X : 0.033,0.02, and 0.01. The remaining input data come from the standard set. 
The dashed lines represent the lcr variation about the current experimental results for the 



individual Z leptonic widths, T e T x T p = 83.85 ± 0.29 MeV and T e e x e p = 83.92 ± 0.17 MeV [gl 
The one-loop SM prediction is r rr = T ee = 84.03 MeV. 

We have discussed the main features of Figs. [7|a,b {Ta ~ rr^Mjf, rising T ee vs falling 
T Tr for TT m i X dominating ee m i X ) previously 0. Here we point out that the main difference 
between Figs. [7|a,b and Figs. |8|a,b, namely the total upward shift of the widths in the latter, 
can indeed be traced to the tree-level correction to the //-decay. Unfortunately, this tree- 
level correction interferes destructively with the one-loop corrections which drive T TT widths 



down. Also note that as TT mix and ee mix become comparable, so do T TT and T ee , as expected. 

Our best constraints at the 2a level on NHL mass come from T TT , shown in Fig. [F]a Q we 
get M N < 4.3 TeV for rr mix = 0.033, M N < 7 TeV for rr mix = 0.02, and M N < 13 TeV for 
TTmix = 0.01. These constraints can be neatly summarized in the following approximation 

033 

M N < 4.3 x TeV. (40) 

TT~ m i x 

Here, this assumes TT mix dominates ee m i X . 

We note that as the value of TT mix is more tightly constrained, these limits are less 
restricted than those from perturbative unitarity considerations pi: 



/0 033 

M N < 4 x TeV. (41) 

V 7~T m i x 

Finally, we note that our computation of the muon decay loops also enables us to find 
the one- loop modification of Eq. [7| by NHL's: 

T{r^evv)/T SM {r^evv) ^2 1 - rr mix 1 + 25? 



T(ji -> euu)/T SM (fi -> evv) ^g^ 1 - [i[i mix 1 + 25^ ' ^ 2 ' ) 
where 5& is given in Eq. For = 4 TeV, and the current constraints on the mixing 
parameters, we find that the one- loop correction is only about 1% of the tree-level correction; 
therefore the constraints of Ref. |y| are indeed independent of the NHL mass. 



IV. FLAVOUR-CONSERVING VS FLAVOUR- VIOLATING DECAYS 

In this section, we review the constraints on the parameters in our model as derived 
from flavour-violating and flavour- conserving decays (the latter will be represented by the 
leptonic decays of the Z boson). We compare the sensitivity of these two classes of processes 
to the presence of NHL's. 

Lepton flavour-violating decays have so far received a lot more attention PJ7|-[14]| than 



the flavour-conserving processes |I|p2]|. The calculation of the flavour- violating processes is 



Constraints from T ee shown in Fig. 0b are just slightly worse at the 2a level. 



simpler, with a smaller number of contributing diagrams and without the need to renor- 
malize. Also there could be a certain preconception that the experimental signature of the 
flavour violation is more 'dramatic'. It is our intention to show here that in many cases this 
expectation is not justified. We give also the results of our calculation of flavour- violating 
decays of the Z boson. A summary of experimental limits, theoretical predictions and the 
constraints on the mixings and/or NHL masses implied by flavour- violating decays is given 
in Table |TJ. We will now address these decays one by one. 

In the case of the flavour-violating mixing parameters, the constraint on one of them, 
e fornix, arises from the measured limit of a rare decay \i — > ej. This decay was studied in the 
context of our model and of see-saw models with enhanced mixings (an example of a see-saw 
model with enhanced mixings is the model of Ref. ||23|| ) by several authors [[7|-[TT|] . The /j, — > 
branching ratio goes like |e/i m j X | 2 times a function which is independent of for > 500 



GeV. The current experimental limit on the \x — > ej branching ratio, BR exp < 4.9 x 10 



-li 



0, yields a very stringent upper limit on the mixing of \e\x m i X \ < 0.00024 (see Table H). 

One might expect the other flavour- violating mixing parameters, er mix and fiT m i X , to 
be limited by the corresponding flavour-violating r decays. However, experimental limits 
on r -> e7 and r — > pq y BR exp < 1.2 x 10~ 4 , BR exp < 4.2 x 10 -6 respectively ||, are 
much weaker. Moreover, the predicted rate BR t h = 7 x 10 -7 || is now out-of-date due to 
improved constraints on the mixings (see Table [II]). With the current limit {rr m i X = 0.033) 
the predicted rate would be smaller by at least one order of magnitude, implying that the 
theoretical result is two orders of magnitude below the experimental upper limit for /ry mode 
and about three orders for mode. As a result, it is not these flavour- violating processes 
which place the strongest limits on the flavour-violating mixing parameters. Rather, for 
fJ'Tmix and er mix , we have to use indirect limits obtained by combining the global analysis 
results for the flavour-conserving mixing parameters with the inequality Eq. [5]. 

Several other flavour- violating processes at very low energies have been considered. An- 
other well-constrained muon decay mode is /x — > e~e~e + (BR exp < 1.0 x 10 -12 , 0), studied 
in Refs. |^^| JiO|JTT| . The calculation shows the quadratic nondecoupling which we will en- 



counter in the lepton flavour-violating decays of the Z boson below. Ref. [11] gives (with 
an assumption discussed therein) the following constraint on NHL mass as a function of 
ee mix , efi mix (see Table 0): 

M% < 0.93 x 1(T 5 r- (43) 

&&mix | &[^mix \ 

Also considered in Refs. |7|,|8|,[TT] is \i — e conversion in nuclei, fi~(A,Z) — > e~(A,Z). The 
resulting constraint on the product ee m i X \efi m i X \ |TTJ is similar to the one above. 

For the flavour- violating decays of the tau into three leptons (r — > e~e~e + , e~fi~fi + , 
etc) we know of no calculation studying the large (TeV) NHL mass limit in the context of 
our model. However, within the see-saw model of Ref. [p3| , Pilaftsis |~4| predicts with the 
current limits on mixings (TT mix = 0.033) and for M N = 3 TeV rates BR th {r — > e~e~e + ) = 
5 x 10~ 7 , BR t h(T — > e~[i~fi + ) = 3 x 10~ 7 which are well below the current experimental 
limit BR exp < 1.4 x 10" 5 (see Table H). 

Finally, hadronic decay modes of the r lepton, r — > Irj, In IJ are also disfavoured by 
loose limits, e.g. BR exp (r ^-yr ) < 4.4 x 10" 5 0. 

Consider now the case of the flavour-violating leptonic decays of the Z boson. These 
rare processes were studied in the context of our model previously P,|T2fl ; however, the limit 
of large NHL mass was not fully investigated. This point was noted in Ref. [13|], where the 
branching ratios for Z —>■ l^l^ (e ± /i =F , yU ± r T , e^-r^) were derived in the see-saw model of 
Ref. [ 23| . We therefore here present the results in our model, having carefully treated the 
case of a large NHL mass but without showing the calculational details. The loop diagrams 
involved are very similar to those of the flavour-conserving leptonic decays of the Z boson 
which we discussed fully in Ref. M and the calculation of the flavour-violating process is 



very similar to that of Ref. [pl| for the other model. The particular predictions depend, 
along with the NHL mass, on the various mixings and their relative phases. We do not 
present full details here since the results are not particularly promising. The branching 
ratio BR(lfl^) = r ; + r+r ^+/r^ is shown as a function of M/v in Fig. |^a,b. The parameter 
5 introduced in these figures corresponds to the allowed range of relative phases arising in 



the mixing factors. In Fig. [|a, we set 5 = — 1 and in Fig. ||b, 5 = +1. 

As an example for comparison with experiment, using maximally allowed mixings 
( TT mix = 0.033) and 5 = +1, we predict the following branching ratio limits for M N = 5 TeV 



(44) 



BRth 


[z - 




BRth 


[z- 




BRthi 


z - 





These results are similar to those of Ref. 14 , where, as noted above, the calculation was 



done in the context of a see-saw model with enhanced mixings. For experimental limits see 
Table f|. Our most promising prediction, for the er mode, is at least one order of magnitude 
below the experimental limit. Hence the flavour- violating leptonic decays of the Z boson do 
not represent a good chance for finding evidence of NHL's. 

We conclude that the flavour- violating processes give only one (mixing dependent) con- 
straint on NHL mass coming from fi — > eee (or \i — e conversion in nuclei), see Eq. [E| For 
ee m j X = 0.0071 and \efi mix \ = 0.00024, this yields M N < 2.3 TeV. For the remaining flavour- 
violating processes to become sensitive to NHL mass, the experimental upper limits would 
have to be pushed down by at least one order of magnitude for flavour-violating leptonic 
decays of the Z boson, and by one to two orders of magnitude for flavour-violating decays of 
the r lepton. This most likely requires increased high luminosity running at LEP I energy 
and a r factory |24| . 



On the other hand, the flavour- conserving processes lead to limits on summarized 
in Eq. HOI, which for TT mix = 0.033 give < 4.3 TeV. These limits depend on different 



mixing parameters than the flavour-violating constraint and thus probe a different part of 
the mixings vs NHL mass parameter space. A disadvantage of the flavour-violating decays is 
that they are always proportional to a flavour- violating parameter and this can lead, via the 



2 The extremely useful limit on e\x mix arising from /i — > ej, is not sensitive to M^, for > 500 
GeV. 



inequality Eq. ^ to their further suppression with respect to flavour- conserving processes. It 
is perfectly possible that there might be signatures of the flavour-conserving processes even 
if there is no sign of the flavour-violating ones. 



V. CONCLUSIONS 

In this paper, we have generalized our previous analysis of a model containing NHL's by 
relaxing the restriction on mixing parameters ee m i X = fifi m ix — 0. This involved evaluating 
one- loop corrections to the muon decay which feed into the input parameter M\y. We found 
that two box diagrams exhibit quadratic nondecoupling but that they are only marginally 
important numerically. Hence the numerical results of Ref. [U remain basically valid, al- 
though they can be improved by the inclusion of the tree-level correction to the muon decay, 

The mass M^, if larger than Mz, can presently mainly be probed in radiative corrections 
(loops). A traditional approach was mostly limited to hypothetical lepton flavour- violating 
processes such as /i — > ej, fi, r — > ee + e~, Z — > e ± ji T 1 etc P,[7|-[14|. We reviewed constraints 
from these processes in Sec IV. 

NHL's could also induce (again via radiative corrections) deviations from the SM in 
currently observed processes, such as those we have previously considered: the leptonic 
widths of the Z boson Tu, lepton universality breaking parameter [/& n and the mass of the 
W boson M\y- The effect of the NHL mass Mjy in such radiative corrections is, on the one 
hand, suppressed by small mixings; on the other hand it is enhanced due to nondecoupling, 



the violation of the Appelquist-Carazzone theorem ||17|| . These competing tendencies are 
reflected by the typical behaviour of the dominant terms, 

{rr mix ) 2 M 2 N . (45) 



To make up for the small mixings, only NHL's with masses in the TeV range can lead to 
significant deviations from the SM. In the case of one mixing, TT m i X , dominating, we found 
(see Eq. |4(|) the following approximate dependence of M N on TT mix {2a level): 



n 033 

M N < 4.3 x TeV (46) 

TT~mix 

which arises from the consideration of Z leptonic decays. We also found some sensitivity 
of the W mass to NHL mass and mixings, but these are quite dependent on the top quark 
mass so we cannot summarize them in the same way. 

These limits on are only matched by those from \i — > eee.The flavour- violating decay 
rates for r, which we reviewed in Sec. IV, and for the Z boson, derived in Sec. IV, are 
below the current experimental sensitivity. Moreover, the /i — > eee decay depends only 
on ee m i X and efi m i X , two of the six mixing parameters, and may be unobservable if ee m i X 
and/or e/i m j X are very small. The inequality Eq. |5| can further suppress the flavour- violating 
processes against the flavour- conserving ones via the 'conspiracy of the phases' in the sum 
of complex terms making up the flavour- violating parameters. 

For these reasons, the first signatures of neutral heavy leptons could come from flavour- 
conserving observables. At this time, LEP has stopped its runs at the Z-peak energy and is 
running at 130 — 140 GeV. It will eventually be producing W pairs which will allow the mass 
M w to be measured with a precision of 0.044 GeV [E§] (currently M w = 80.410 ± 0.180 



26f). Combined with more precise measurements of the top quark mass we might be in a 
position to place even more stringent limits on NHL masses and mixings from our prediction 
ofM w . 
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APPENDIX: A 



The total contribution of the box diagrams (Figs. [I] a-g) is 



■M-box = -M-ZeWfi + ■M-WuZu + -M-WvZN + M-WNZv + M-WNZN 
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The computation of the vertex diagrams (Figs. 0a-f) yields 



Vertex = M^Z + M^z + M 
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The part of the charged lepton self-energy which we specifically use in the text is: 
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The left-handed part of the neutrino self-energy (Fig. |5|) is given by 
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APPENDIX: B 



Here we prove Eq. [22], A^ HN + A^ xN + 5Zf N = 0, using Eq. |21j 



We note the vertex V^ N (Fig. |c) is given as 



(2n)"s/2M w y 2 h- A-M, 



N 



V2M W K H)aX N 2 ( q - Pl - P2 )2- M ^q^-M^ 
x i e (-2g + p 1 +p 2 ) (Bl) 

If we now return from the 7// vertex to the W fxv vertex (Fig. §f), we have 

yr W v f d n q -ig 2 f ! + 7s j 

*W ( 27 r)«2M lv Mjv( ^^ )ia 2 A-fr-M N 

V2M W ^ K ^ Mn 2 ( g - Pl - p 2 y -Ml?- Ml 
x ^( + 2q - Pl -p 2 ) . (B2) 

A similar expression holds for VY^ N . The Lorentz structure of Eqs. [Bl] and [B!| is the same. 
In the large M^ limit, M# and My/ in the propagators are negligible; therefore the only 
possible difference between the two vertices comes from constant factors. If we forget for a 
moment about the mixing factors, it can be easily checked that 

v w _ Vjjgg vl 

v <pHN — ^ v 4>(j>Ni 
M.<j>HN + -MfaN = (A^HN + ^4>xN)-M-tree 

A^jv + A^* = 2Ff w = 2Ff JV . (B3) 
Hence (using Eq. [B3] and Eq. |2l]) 



= 2(^^ + 5^) = 0. (B4) 

That is, the two dominant nonstandard contributions from Eq. [16] cancel. To show that the 
inclusion of the mixing factors will not affect Eq. IB3, note that the mixing factor for the 



W [iv vertex, which we denote as k±, is related to that of the 7// vertex, denoted as k 2 , in the 
approximation in which flavour- violating mixing factors e/i mix , Tfi mix , er mix are vanishing, as 
follows. 

a 

= (41(^(4),, + (^E(^A + (^irE^U^ 

a a a 

= {K[) ie enmix + {Kl)^fifi mix + {Kl) iT Tfi mix 

= {K\\^ mix = (K[)^k 2 (B5) 
The remaining factor, (K^)., is absorbed into A4 tree as required. 
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FIGURES 



FIG. 1. Box diagrams for muon decay. 
FIG. 2. Vertex diagrams for muon decay. 
FIG. 3. Charged lepton self-energies. 
FIG. 4. 7// vertex . 
FIG. 5. Neutrino self-energy diagrams for muon decay. 

FIG. 6. Counterterm diagram for neutrino self-energy in muon decay. 

FIG. 7. Z leptonic width as a function of Mn for ee m i x = 0, m t = 176 GeV, Higgs mass = 200 
GeV and different values of the mixing parameter TT m i x , (a) Z — ► rr mode, (b) Z — > ee mode. The 
dashed lines represent the la band about the current experimental value (a)T^ p = 83.85 ± 0.29 
MeV, (b)r^ = 83.92 ± 0.17 MeV. 

FIG. 8. Z leptonic width as a function of Mn for ee m i x = 0.0071, m t = 176 GeV, Higgs 
mass = 200 GeV and different values of the mixing parameter TT m i x , (a) Z — ► rr mode, (b) 
Z — ► ee mode. The dashed lines represent the 1<t band about the current experimental value 
(a)r^ = 83.85 ± 0.29 MeV, (b)r^? = 83.92 ± 0.17 MeV. 



FIG. 9. The branching ratio Z — > ifl^ as a function of Mn for (a) <5 = — 1, (b) 8 = +1. 



TABLES 





SM 


M N = 0.5 TeV 


M N = 5 TeV 


M N = 15 TeV 


M N = 30 TeV 




ee mix = 0.0071, /ifj ma = 0.0014, TT mix = 0.0 
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ee mix = 0.0071 


, WW = 0.0014, TT mix = 0.033 
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TABLE I. Contribution of the muon decay loops to <5y, A r and Mw- 



Process 


Experimental BR H 


Theoretical BR 


Limits on masses and/or mixings 


H -> e7 


< 4.9 x KT 11 90% C.L. 


4.9 x 10~ n a 


|e/w| < 0.00024 


r — > e7 


< 1.2 x 10" 4 90% C.L. 


7 x 10" 7 b 
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r -> ^7 


< 4.2 x 10~ 6 90% C.L. 


7 x 10~ 7 b 
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/u — > eee 
r — > eee 


< 1.0 x 10~ 12 90% C.L. 

< 1.4 x 10~ 5 90% C.L. 


1.0 x 10~ 12 c 
5 x 10~ 7 d 


Mjy < 0.93 x 10~ 5 1T , ey2 | 




r — > eii/x 


< 1.4 x 10~ 5 90% C.L. 


3 x 10~ 7 d 




Z — > e/x 


< 6.0 x 10~ 6 95% C.L. 


3.3 x 10~ 8 e 




Z — > er 


< 1.3 x 10~ 5 95% C.L. 


1.4 x 10~ 6 e 




Z — > /iT 
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2.2 x 10~ 7 e 
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| = 0.00024, M N > 0.5 TeV; Ref. § 








= 0.043, = 0.008, TT mix = 0.1, M N > 0.5 TeV; Ref. g 




c M 2 N x 


ee m «|e/x mi:c | = 0.93 x 10~ 5 x 1 TeV 2 ; Ref. 







= 0.01, fJ,n m ix = 0, rr mil = 0.033, M N = 3 TeV 


; Ref. @ 






= 0.0071, /i/w = 0.0014, rr m!I = 0.033, M^y = 


5 TeV; this paper and Ref. 
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TABLE II. Flavour-violating decays: experimental limits, theoretical predictions and the con- 
straints implied. 
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